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pH . Abstract 
• 1—1 ' 

I We compute accurately the golden critical invariant circles of several area-preserving 

I— I' twist maps of the cylinder. We define some functions related to the invariant circle and 

I . to the dynamics of the map restricted to the circle (for example, the conjugacy between 

^ I the circle map giving the dynamics on the invariant circle and a rigid rotation on the 

■ circle) . The global Holder regularities of these functions are low (some of them are not 
. even once differentiable) . We present several conjectures about the universality of the 

0> I regularity properties of the critical circles and the related functions. Using a Fourier 

^ ' analysis method developed by R. de la Llave and one of the authors, we compute 

. numerically the Holder regularities of these functions. Our computations show that 

i - withing their numerical accuracy - these regularities are the same for the different 

I maps studied. We discuss how our findings are related to some previous results: (a) 

^ ■ to the constants giving the scaling behavior of the iterates on the critical invariant 

^ , circle (discovered by Kadanoff and Shenker); (b) to some characteristics of the singular 

'k^ I invariant measures connected with the distribution of iterates. Some of the functions 

■ studied have pointwise Holder regularity that is different at different points. Our results 



give a convincing numerical support to the fact that the points with different Holder 
exponents of these functions are interspersed in the same way for different maps, which 
is a strong indication that the underlying twist maps belong to the same universality 
class. 
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Area-preserving twist maps of the cylinder are popular models of the dy- 
namics of many physical systems, i.e., they occur as Poincare maps of 2-degree- 
of-freedom Hamiltonian systems. Homotopically non-trivial invariant circles of 
such a map play an important role in organizing the global dynamics of the 
map. Generally, as the perturbation grows, more and more of these circles are 
destroyed, until there remains only one such circle, called the "critical" circle. 
This circle is the last obstacle to the unbounded growth of the "action variable". 
In this critical situation many characteristics of the system become drastically 
different from the "under-critical" case. For example, consider the dynamics of 
the iterates of the twist map restricted to the critical circle — it is given by a 
map of the circle. This map can be conjugated to a rigid rotation on the circle, 
but the conjugating function has very low regularity — its Holder exponent is 
lower than 1. The Holder regularity of this conjugacy is related to some uni- 
versal properties of the map, i.e., to the universal rescaling factors fP, '2j and 
to the scaling properties [3J of the distribution of the iterates on the critical 
circle (which is governed by a singular invariant measure). We compute several 
functions related to the critical circle and to the dynamics on it and us a method 
developed in [4j to assess numerically their global Holder regularity. Our find- 
ings lend support to several conjectures concerned with the universality and the 
renormalization group description of these critical objects. 



1 Introduction 

It has been known since the late 1970's and early 1980's that many objects at the boundary of 
chaotic behavior exhibit remarkable scaling properties and that, furthermore, these proper- 
ties are "universal". Such properties are exhibited by unimodal maps of the interval 00 EI, 
critical maps of the circle [H], critical KAM tori 0121, and other systems. These observations 
were explained in terms of a renormalization group analysis, following a methodology that 
had been developed earlier in the study of critical phenomena in statistical mechanics and 
field theory [TOl HH [12] . 

The scale invariance of the critical objects affects many of their properties. Notably, the 
Holder regularity of the critical objects (or some functions related to them) tends to have a 
low and fractional value. Presumably the values of the regularities are related to exponents 
and geometric properties of the renormalization group fixed points which describe the critical 
objects. 

Furthermore, the observation that critical objects can be divided in "universality classes" 
such that all objects in a given class "look the same" can be tested numerically. One way to 
do this is to define certain functions related to the critical objects - typically these functions 
are not very regular (in some cases not even once differentiable), and to test numerically 
whether the regularities of these functions are the same for different objects. Another - even 
more sensitive - test for universality is to take two functions, say hi and h2, from the same 
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class, and to study the regularity of functions like hi o - for hi and /i2 belonging to 
the same universality class, one can expect that hi o Z?,"^ be more regular than hi and ^• 
These ideas were tested in ^ in the case of non-critical and critical (with different degree 
of criticality) circle maps, in which the empirical results are accompanied by an extensive 
mathematical theory. A substantial part of the effort in |3] was to develop implementations 
of methods known in harmonic analysis (finite differences, Littlewood-Paley theory, wavelet 
analysis) to assess the regularity of the objects numerically. 

In the present paper, we extend the methodology of |3] to the study of critical invariant 
circles of area-preserving twist maps. Invariant circles in dynamical systems are among 
the most important objects that organize the long-term behavior of the system, and the 
critical ones are especially important because of their role as "last barriers to chaos" (for 
readable reviews see, e.g., ^M,-, or, with more emphasis on the mathematical aspects, the 
recent book Critical invariant circles have been extensively studied since early 1980's 

P 121131121 • 

We compute accurately the golden critical invariant circles of several standard-like area- 
preserving twist maps and some functions related to the dynamics of the iterates of the maps 
on these circles. Then we apply methods developed in jl] to study the Holder regularity of 
these functions and some universality aspects. 

In Sec. m we give some background on twist maps and their critical invariant circles, 
define the functions that are the objects of our numerical study, and present several precise 
conjectures concerning the properties of the critical invariant circles and the functions intro- 
duced. Sec. El is devoted to a discussion of the numerical methods used to compute critical 
invariant circles and to assess Holder regularity of functions. We collect our results in Sec. EJ 
and in Sec. Owe discuss their significance and relationship with previous studies. 

2 Critical invariant circles of twist maps 
2.1 Twist maps 

Let T := M/Z stand for the circle. We will be concerned with maps F of the (infinite) 
cylinder T x M, 

F : T X M ^ T X M : (^, r) ^ F(0, r) =: {6', r') , 
which satisfy the following properties: 

• Area preservation: The map F preserves the oriented area: deiDF = 1. 

• Zero-flux: The oriented area between a homotopically non-trivial circle and its image 
under F is 0. (In our situation, this is equivalent to saying that every non-trivial circle 
intersects its image.) 

• Twist condition: For any fixed value of 6', ^ > 0. 
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A map of the cylinder can be identified with a map F : ^ (called a lift of F) 
which satisfies 

F{e + 1, r) = F{e,r) + (1,0) . 

Often one does not need to keep the distinction. 

The maps which we will use in our numerical studies are of the form {6',r') := F{6,r) 
with 

e' = {e + r') mod 1 , 

(1) 

r' = r + \V{e) , 

where A is a parameter, and : T ^ R is a function satisfying V{6) dO = 0. In particular, 
many numerical studies have been devoted to studying (0) with 

V{e) = — ^sin27r^, (2) 

271 

in which case we will call the map F the Taylor- Chirikov map. 

Given an orbit X = {{On, f^n) = -^"(^O; ''"o) I = 0, 1, 2, . . .}, we define its rotation number, 
p{X), as the limit 

p{X):= hm 



*±oo n 

whenever this limit exists. In contrast with the situation for circle maps, the rotation number 
depends on the orbit (and it may happen that some orbits do not have rotation number). 

We say that an orbit is well- ordered when for every k and /, the function of n defined as 
e(n) = 6n+k — l — On has the same sign. Every well-ordered orbit has a rotation number (the 
converse, however, is not true). 

It is also easy to see that if a bounded orbit is well-ordered and p(6'o,ro) is irrational, 
the closure of the orbit, {(6'„, r„)}5^Q, is a perfect set (i.e., every point is an accumulation 
point of points in the set); in other words, in this case the orbit is either a homotopically 
non-trivial circle or a Cantor set. 

A set t/ C T X R is invariant liU = F{U). 

The following result plays an important role. 

Theorem 2.1. // F is as above, for every p G R there exists a well-ordered orbit with 
rotation number p. 



2.2 Invariant circles of twist maps — rigorous results 

The proof of the following theorem can be found in , ^Z] • We refer to and jTHj for a 
detailed exposition. 

Theorem 2.2. Let U be an open simply connected invariant set, containing one of the ends 
of the cylinder. Then the boundary, dU , of the set U is an invariant circle which is the graph 
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of a Lipschitz function. In other words, dU can be written as r = R{9), where i? : T — R 
a Lipschitz function. 

For the map (Q), the Lipschitz constant of the function R can he hounded hy an expres- 
sion which involves only the Lipschitz constant of the function F in a neighhorhood of the 
circle dU . 

In particular, 

Corollary 2.1. Any homotopically non-trivial invariant circle is the graph of a Lipschitz 
function R. In the particular case of the map ((T}, the Lipschitz constant of R can he hounded 
hy a constant which is independent of XLipV. 

A number p is said to be Diophantine if, for each m,n G N \ {0}, for some C > 0, and 
for some d > 2At satisfies 



m 

P 

n 



c 

r? 



In the case when the map F is close to integrable and its rotation number p is Diophantine, 
one can apply Kolmogorov- Arnold- Moser theory to obtain that there exists an analytic 
invariant circle such that the orbits on it have rotation number p. 

Golden invariant circles are those with rotation number equal to the golden mean, 

crG:= [1,1,1,...] = ^^. (3) 

Here we have used the notation p = [ai, 02, 03, . . .] = l/(ai + l/(a2 + l/('^3 + ■ ■ ■ ))) ^^e 
continued fraction expansion of p G (0, 1) [19 . 

There are also rigorous results that guarantee the non-existence of invariant circles of F 
of the form ((H). 

Theorem 2.3. (i) //supg |AV^(^)| > 1, then has no invariant circles. 

(a) If sup \ V'{9) \ = 1 (which holds for the function then for |A| > | the map ((H) has 
no invariant circles. 
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(Hi) For V given hy the map (0) has no golden invariant circles for |A| > 
0.984375.... 

(iv) For V given hy (j2]), the map has no golden invariant circles for |A| > 0.9718. 

Part (i) of Theorem 12.31 is elementary: if A supg |^(^)| > 1, then there will exist points 
(r,r*) G T X M such that F{e*,r*) = (r,r* + 1), which, iterated, gives that F"(r,r*) = 
{6*,r* + ra) - the unbounded growth of the second coordinate of F"'{6*,r*) with n implies 
that a topologically non-trivial invariant circle cannot exist. 

Part (ii) can be found in , parts (iii) and (iv) are proved by computer-assisted methods 
in |2ni and EH, resp. 

It is widely believed that 
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Conjecture 2.1. For a Diophantine number p and for a map F of the form (Q, there is a 
number A{p) such that when |A| > A(p), there is no invariant circle with rotation number p, 
and when |A| < A{p), there exists an analytic invariant circle with rotation number p. The 
invariant circle becomes critical when |A| = A(p). 

Since our paper is devoted to homotopically non-trivial invariant circles, we will usually 
omit the words "homotopically non-trivial" . 

2.3 Functions related to the critical invariant circles 

We are interested in describing the critical invariant circles with rotation number p which 
are in the boundary of existence. Postponing for the moment issues on how these objects 
can be actually computed, we point out that to a given critical invariant circle 7 of rotation 
number p, we can associate: 

• the function i? : T ^ M such that the critical invariant circle 7 is the graph of R: 

7 = {(^, r) G T X M : r = R{e)} ; (4) 

• the advance map g : T ^ T defined by 

F{e,R{e)) = {g{e),Rog{e)) (5) 

• the hull map : T ^ T x M, which gives a representation of the invariant circle 7 in 
such a way that the dynamics on 7 becomes a rotation by p, i.e., 

F o^{e) = "^{9 + p) ; (6) 

• the map /i = 7rio\l/:T^T (where tti : T x R ^ T is the projection onto T), which 
conjugates the advance map to a rotation by p: 

g o h{e) = h{e + p) ■ (7) 

• the map : T ^ T, which is the inverse of the map h defined in 
We note the following rigorous results. 

Theorem 12.21 guarantees that the invariant circle the function R is Lipschitz. It is an easy 
consequence of the implicit function theorem that g should be as regular as R. Nevertheless, 
it is useful to compute the regularities of both g and R independently to asses the reliability 
of the numerical methods used. 

Because of © , it is clear that the regularity of g is not smaller than the minimum of the 
regularities of h and h~^. 



6 



2.4 The "big" conjugacies 



Let p be a Diophantine number, Fi {i = 1, 2) be area-preserving twist maps, and 7, be 
the critical invariant circle of Fi with rotation number p. Let g^^ and h.y. be the associated 
advance map (0) and conjugacy resp. We introduce the conjugating functions 

G'71,72 := fi'71 oS-^a^ : T ^ T , 

We will call these functions "big" conjugacies to distinguish them from the "small" conju- 
gacies h that conjugate the projected dynamics on the critical circles to a rigid rotation (jZj). 
Note that the "big" conjugacies satisfy 

^71,72 ° ^72,73 ^71,73 ' -^71,72 ° -^72,73 -^71,73 • 

Below we discuss one aspect of the definition of the big conjugacies that will be important 
in our computations. 

Since there is no "origin" on the circle T, one has certain amount of freedom in the 
definition of some maps. For example, if the function \E' is a hull map (i.e., satisfies (jH))), 
then the function defined as ^(^) = "^{0 + () will also satisfy (jH)) for any choice of the 
constant (. Similarly, the map h ((Zj) that conjugates the advance map g to a rigid rotation 
can be redefined by composing it on the right with a rotation, and the resulting map, 
h{6) = h{6 + (), will also conjugate g to a rigid rotation. Naturally, all important properties 
of the maps h and h - in particular, their Holder regularity - will be the same. However, 
one cannot use this freedom liberally when studying the big conjugacies. To understand the 
reason for this, consider the map h defined by ((Tj) for some twist map F. Naturally, the 
map h o is the identity map, so it is C°°. However, for any nonzero ( in the definition of 
h, there is no guarantee that the map h o h^^ will be C°°. This is due to the fact that the 
regularity of h may be different at different points, and while in hoh~^ these "irregularities" 
cancel out, in h o the action of h does not necessarily "undo" the irregularities caused 
by . In Sec. 12.51 we explain in detail how we choose C, in order to avoid the "spurious" 
irregularities of the big conjugacy. 



2.5 Big conjugacies and symmetries 

Consider two functions h^^ and h^^ corresponding to the critical circles 71 and 72 of the 
twist maps Fi and F2. If Fi and F2 happen to belong to the same "universality class" (see 
Sec. 12. 6p . then one would expect that the big conjugacy H^^^^^ will be more regular than 
the functions h^^ and h~^. To avoid introducing spurious irregularities in H^-^^^^, we use the 
symmetries of the map h that come from the symmetries of the F |^ • 

It is well known that if the function V is odd, then the map F given by can be written 
as a composition of two involutions: 

F = Ji o Jo , J2 = /2 = Id , (8) 
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where 

Io{e,r) = {-e,r + \V{e)) , h{e,r) = {-e + r,r) . (9) 

From (jHl) we have Iqo F = F^^ o Jg and Ii o F = F^^ o Ii. Acting on (jHl) with Iq from the 
left, we obtain 

F-io(/oo^)(0) = (/oo^)(0 + p) . (10) 

On the other hand, if we define the function L : T ^ T x M by L{6) := ^(—0), then © can 
be written as 

F-' o L{9) = L{9 + p) . (11) 

Comparing (fTUI) and ()11|) . we see that L and Jq o \[/ can differ only by a shift in the argument, 
i.e., there has to exist a constant ( such that Iq o \E'(6') = L{6 + C) = "^{—0 — (). This, 
together with Q and /i = tti o \[', imply 

/i(^) = _/i(_^ _() . 

This implies that h{—^) = 0, and the numerical value of ( can be found from the computed 
values of h. Setting h{6) := h{6 — |), we obtain that h is an odd function. In what follows, 
we will assume that the appropriate value of ( has been subtracted, and will omit the tilde 
over h. 

2.6 Universality 

In this section, we formulate precisely some conjectures on the behavior of critical invariant 
circles described by a non-trivial fixed point of the renormalization group. It seems quite 
possible that these conjectures can be proved as conditional theorems assuming existence 
and certain properties of this fixed point. 

One of the most striking predictions of the renormalization group theory is that many 
characteristics of the critical invariant circles are largely independent of the details of the 
map. This is captured by the notion of universality. 

Definition 2.1. We say that a numerical characteristic is universal when it takes the same 
value in an open set of functions. We say that a property is universal when it holds for an 
open set of functions. 

The open sets alluded to in Definition 12.11 are called domains of universality. 

For the case that we will be concerned with, the description of the domains of universality 
in terms of properties of the non-trivial fixed points of the renormalization operator is still 
debated, but there are indications that the domain of universality is not the whole space 

[2nii2niizzi- 

Conjecture 2.2. The existence of one and only one non-trivial fixed point of the renormal- 
ization operator is a universal property. 
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This conjecture has been known for a long time jT^j. Recently in |2H1 it has been shown 
how this conjecture follows rigorously from an extension of the standard renormalization 
group picture. Even the formulation of the subsequent conjectures depends of Coniecture l2.2l 

The concept of universality is rather natural when one wants to study properties that 
depend on the speed at which the set of maps converges to the fixed point under the renor- 
malization operator. In particular, regularity of conjugacies depends on the this speed of 
convergence and, hence, should be a universal quantity (more precise formulations are given 
in [211 )■ Hence, we can also conjecture that: 

Conjecture 2.3. The regularity, n{R), of the critical invariant circle is a universal number. 

Conjecture 2.4. The regularities n{g), K,{h) and K,{h^^) are universal numbers. 

Conjecture 2.5. For pairs of critical circles and'~f2, the regularities n{G^f^^^^) and nlH^-^^y.^) 
are universal numbers. 

Directly from the definition of Holder regularity, one can see that if K{(f)) and are 
between and 1, then K{(f) o ip) > K,{(f)) This implies that 

'^(^71,72) = i^iK ° K^) > i^iK^) ■ (12) 

For all critical invariant circles 7^ that we studied, we obtained numerically that n{h^J < 1 
and i^{h~l) < 1, so (fT^ yields that H^^^y^ is not less regular than i^{h^^) K{h~^). For 71 and 
72 in the same universality class, however, we expect more - because of "cancellation" of the 
"singularities" of h^^ and , we state our final 

Conjecture 2.6. The following inequalities hold for « = 1, 2: 

l^iK) < '^(^71 ,72) . '^(^7/) < '^(^71 ,72) • 

3 Description of the numerical methods 

In this section we first describe the methods used for numerical computation of invariant 
circles and the related functions described in Sections 12.31 and 12.41 Then we briefly discuss 
the method we use to compute the global Holder regularity of the functions. 

3.1 Computing critical invariant circles 

We need to compute (homotopically non-trivial) critical invariant circles of twist maps of 
the form with a Diophantine rotation number. We approximate such invariant circles 
by well-ordered periodic orbits (whose existence is guaranteed by the so-called Birkhoff's 
Geometric Theorem [20] )■ Consider a sequence {X^^^j^^q of well-ordered periodic orbits 
whose rotation numbers, {pj}jm, constitute a sequence of rational numbers which converge 
to a Diophantine number p. Then the limit of these periodic orbits will be a well-ordered 
invariant set Xp of rotation number p; the existence of this set is guaranteed by Aubry-Mather 
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theory Ch. 13], [TU Ch. 2]. The set Xp can be a continuous curve which is a graph of 
a Lipschitz function under appropriate conditions (Theorem 12. 2p or an orbit homeomorphic 
to a Cantor set (Cantorus). 

We approximate a Diophantine number p by the rational numbers given by finite trunca- 
tions of the continued fraction expansion of p. In the case of the golden mean ctq these 
rational approximants are ratios pm = Qm-i/Qm of consecutive Fibonacci numbers Qm- The 
limit of the periodic orbits with rotation numbers pm is the invariant set Xp we are looking 
for 123 • 

The problem of computing well-ordered orbits with a prescribed rational rotation number 
Pm is greatly simplified if the function V{6) in is odd. In this case the task of finding a 
periodic orbit is reduced to a one-dimensional problem because the map F can be written as 
the composition of two involutions as in (jHl); if such a decomposition is possible, the map F 
is said to be reversible. If F is reversible, there exists a set of straight lines in the {6, r) space 
- called symmetry lines - that are invariant with respect to the maps Iq and Ii. It can be 
shown that any periodic orbit has two points that belong to one of these invariant straight 
lines, hence we can find these points (and, therefore, the periodic orbits that contain them) 
by using a one dimensional root finder [23j. Using the fact that the periodic orbits computed 
in this way are well-ordered, we can implement a numerical procedure to compute periodic 
orbits of several million points that approximate the invariant set Xp. 

We are interested in studying the properties of area-preserving twist maps of the form 
(^. When the parameter A in is equal to 0, the corresponding twist map acts on each 
point {6,r) as a rigid rotation in 6'-direction, F[6,r) = [6 + r,r), hence the phase space is 
foliated by invariant circles of the form {r = const}. For small values of |A|, KAM theory 
guarantees the existence of invariant circles with Diophantine rotation numbers. According 
to Conjecture 12.11 there is an upper bound A(p) on the values of |A| such that for |A| < A(p) 
there exists an invariant circle with rotation number p; (some rigorous upper bounds on A(p) 
are given in Theorem 12. 3|) . To find an accurate numerical approximation of the critical value, 
A(p), of A for which the invariant circle of rotation number p disintegrates, we applied an 
empirical method known as the "residues method" proposed in [221; developed in [32], and 
partially justified rigorously in The main idea of this method is to determine the value 
of A such that the residue of all the approximating periodic orbits reaches the same value. 
Let Rm be the residue of a periodic orbit which is the mth approximant to an invariant 
circle with rotation number p. If limm^oo Rm = 0, then there exists an invariant circle with 
rotation number p; if lim^^oo Rm = oo, then the invariant set Xp becomes a Cantor set. A 
golden critical invariant circle is obtained at the value of A for which Rm — —2.55426 for all 
values of m. 

3.2 Studying global Holder regularity numerically 

In this section we describe briefly the method we employed to study global Holder regularity, 
referring the reader to ^ for details, references, and assessment of the numerical accuracy 
of various numerical methods for computing regularity. 

In this paper, we will only use the method developed in |^ that was found to be the most 
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accurate for studying global Holder regularity - the so-called " Continuous Littlewood-Paley" 
(CLP) method. Here we do not use the wavelet-based methods implemented in 0. The 
CLP method has been used in I24j . 

3.2.1 Theoretical basis of the CLP method 

We recall the following definition. 

Definition 3.1. For n = n + x with n G Z, % G (0, 1), we say that the function K : T ^ 
has (global) Holder exponent k and write K G Ak(T) when K is n times differentiable and, 
for some constant C > 0, 

\D''K{e) - D''K{e) \ < c\e - e\>' 

for all e,e eJ. 

For the case of an integer value of k, this definition is more complicated, but we will omit 
it since in the applications considered in this paper k is not an integer. 
The following result can be found in Ch. 5, Lemma 5]. 

Theorem 3.1 (CLP). The function K G L°°(T) is in Ak(T) if and only if for some integer 
7] > K there exists a constant C > such that for any t > 

< Ct"-" , (13) 

L°°(T) 

where A is the one- dimensional Laplace operator: AK{6) = K"{6). 

Remark 3.1. If the above result holds for some integer rj > a, then it holds for all integers 
fj > a. 

Remark 3.2. The operator e~*^^~^ is a convolution with the Poisson kernel: e~*^^~^K = 
i^cxp(-2^t) * K. The function u{e,t) := e-^'^K{e) is a solution of Laplace's equation, 
uee + Utt = 0, on the half-cylinder (6*, t) G T x (0, oo), with Dirichlet boundary condition 
u{e, 0) = K{e). 

Remark 3.3. The mathematical theory only requires that (fT^ be an upper bound. In our 
numerical experiments, however, this bound is saturated for a significant range of values of t. 
This fact is very possibly a consequence of the self-similarity at small scales of the functions 
we consider (which is at the basis of the renormalization group description). This saturation 
was also observed for the functions considered in |21] . 
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3.2.2 Remarks on the numerical implementation 

To use the CLP method, we need to apply repeatedly Fast Fourier Transform (FFT), which 
is easiest to do if the values of the function K in ()13|) are known at 2^ equally spaced points 
in the interval [0, 1) for some positive integer N . However, as we describe in Sec. |31 we do 
not have control over the set of points at which the values of K can be computed (where K 
stands for any of the functions i?, h, h~^, H, G). Hence, the first step in applying the 
CLP method would be the computation of the values of K on an evenly spaced grid. If we 
know accurately the values of at M points in [0, 1), we can expect that by using some 
interpolation method, we will be able to obtain the approximate values of on 2^ ^ M 
equidistant points, {2~^j}'^^Q^. To compute the approximate values of K on the equidistant 
grid, we used cubic spline interpolation. Using interpolation poses the question of whether 
the interpolated values represent faithfully the true values of K. Naturally, the answer to 
this question is no, but practically if M is large enough, the interpolated values will be 
very close to the true values, which will allow us to compute many Fourier coefficients of K 
accurately. The degree of "contamination" of the Fourier spectra due to the interpolation 
depends on the uniformity of the distribution of the M points at which the values of K is 
known accurately (see Remark 14. 2p . 

To apply the CLP method numerically, we observe that the operator used 

in Theorem Kill is diagonal in a Fourier series representation: if K{9) = Ylkez^k^''^'"^'''^ , 
then 

'K{e) = ^(-27r|fc|)''e-2'^*l'=li^fee-™. (14) 

fcez 

Having computed the values of the spline interpolant to the function K on an equally spaced 
grid, applying ()13p is easy. Namely, we fix some values of the parameters t] and t, perform 
FFT to find Kk and compute the Fourier coefficients of (^)'' e~*^'^~^i^. Then we apply in- 
verse FFT to find the values of (^-^Y e'^^^'^K at the equally spaced set of points {2^^ j}'^'lo^] 
among these values we find the one with maximum absolute value - this value we take for the 
numerical value of the left-hand side of (fT^ . For a fixed value of rj, we repeat this procedure 
for many values of t (we used several hundred values of t in our computations). According 
to (0, if we plot 

versus logt , (15) 

L°°{T) 

the points should lie below a straight line of slope k — i]. As pointed out in Remark 13.31 fsee 
also Remark I4.3|l the points on the log-log plot should not only be below this straight line, 
but should close to it. We perform linear regression to find the slope of this line, from which 
we find K. 
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4 Numerical results 



4.1 Twist maps studied 

We study numerically a set of one-parameter families of area-preserving twist maps of the 
form each family having a different function V. Within each family we find numerically 
the value A(ctg) of the parameter A for which the golden invariant circle is critical. The set 
of functions V (all of them odd) that we selected is the following: 

1. The standard (Taylor-Chirikov) map: 

V^{0) = -—sm27Te . (16) 

2. The "standard map with two harmonics": 

V2{9) = [sin(27r^) - 0.03sin(67r^)] . (17) 

ZTT 



sin(27r^) - - sin(67r^) 
2 



3. The "critical standard map with two harmonics": 

(18) 

For this choice of coefficients, the first three derivatives of V{6) at 6 = are zero. 

4. The "0.2-analytic map": 

vm = tf''l„, ■ (19) 

2,1x1 — 0.2 cos(27rt^) 

This map has infinitely many nonzero Fourier coefficients. It would be very interesting 
to study this map when the coefficient of the cosine function in the denominator is 
close to 1, but then it would be extremely difficult to compute periodic orbits. 

5. The "0.4-analytic map": 

^'^^^ - "2^1-0.4cos(27r^^) " ^^^^ 

6. The "tent map": 

17 

v,{e) = Y,c,^H'^T^j9) , (21) 

where Cj = (—1)^:;^ for j odd, and Cj = for j even, are the Fourier coefficients of 
the function 

( -49 , for < < i , 
V{9) = \ A9-2, fori<0<|, 
( A-A9 , for I < < 1 . 

The function Vq is close to the piecewise linear continuous function V. 
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Our numerical experiments were performed with the twist maps coming from the above 
six functions V{9) and the corresponding values A((Tg), following the steps below. 

1. As discussed in Sec. Kllf the invariant circle of rotation number ctq can be obtained as 
a limit of periodic orbits of rotation numbers equal to ratios of consecutive Fibonacci 
numbers, Pm = Qm-i/Qm- We chose to compute hyperbolic periodic orbits, and found 
the values of A{cxq) by applying Greene's residue criterion. 

2. The highest approximant to the critical invariant circle that we computed was a peri- 
odic orbit with rotation number Q29/Q30 = 832040/1346269. The value of A(crG) was 
determined by using the condition that the difference, |i?3o — -R29I of the residues of 
the periodic orbits with periods Q29 and Q30 be zero (in practice, we wanted that this 
difference be smaller than 10~^°). The periodic orbits were computed with an error 
not exceeding 10^^^. 

3. We computed the hyperbohc periodic orbit {{Om,rm)}m=o ^f period M = Qso- The 
values of the advance map g 1^ at the points 9^ ijn = 0, 1, . . . , M — 1) were then 
computed by g{9m) = Om+i (here and below, we take modi wherever needed). The 
values of the conjugacy h at the points mac (which corresponds to m applications of the 
rigid rotation by ctq to 0) are given by himao) = 9m, and, similarly, h~^{9m) = maQ. 

4. In our Fourier-analysis-based CLP method we need to deal with periodic functions, so 
we compute the "periodized" versions, g — Id, h — Id, — Id, of the functions g, h, 
and h~^. Then we sort the periodized functions with respect to their argument; the 
function R is already periodic, so we just sort its values. 

5. The periodic functions are passed to the cubic spline interpolation routine to find 
approximations to the values of the corresponding functions on a uniformly spaced 
grid of 2^ points; we used = 20. 

6. The interpolated values of the functions are given to the CLP algorithm to compute 
their Holder regularity. We used integer values of 77 in (jl3|) from 1 to 5, and for each 
analyzed function chose the value of rj that gave the best straight line on the log-log 
plot p5|) . The log- log plots for the other values of rj were used as a consistency check. 

Remark 4.1. In computing the big conjugacies H^-^^^^, we had to take special care of the 
preserving the symmetries of the maps h. For each critical circle we studied, we needed to 
find the appropriate value of the constant ( and shift the argument of the corresponding 
function h as explained in Sec. 12.51 

4.2 Critical invariant circles — visual explorations 

In Fig. n we show the critical invariant circles which, by the definition (0]), are graphs of 
the functions R corresponding to the six twist maps studied. The graphs of the "periodized 
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versions" of the advanced maps, g — Id, the conjugacies, h — Id, and their inverses, — Id, 
are plotted in Figures |21 El andEJ resp. 

Figure illustrates the self-similar nature the functions h; needless to say, the insets are 
true zooms of parts of the graph of the function. 

Fig. El shows the graphs of several periodized big conjugacies H — Id; it is obvious that 
these functions are smoother than the "small" conjugacies h. 

4.3 Fourier spectra, CLP method 

Fig.[7|depicts log^^Q of the modulus of the kth Fourier coefficient of a periodized conjugacy {h— 
Id) versus log^g here h is the conjugacy corresponding to the twist map F with V3 (fTH)) . The 
horizontal distance between two adjacent high peaks is approximately equal to | log^o'^cl ~ 
0.209, which is a manifestation of the self-similarity at small scales. The log^o'l^Sio plots of 
the Fourier spectra of the functions {g — Id) and {h~^ — Id) for the same map F are given in 
Fig. El 

Remark 4.2. Note that the spectrum of h is very accurate even at length-scales ~ 10~^, 
while the spectrum of is quite noisy. As explained in Sections 13.2.21 and 14.11 the main 
reason for this is that the exact values of h are known at the points (mac) modi, which 
are almost uniformly distributed on T. On the other hand, we know the exact values of 
of g and at the very nonuniformly distributed points of the form g"^{6o) (because the 
underlying invariant measure is singular - see Sec. 15. 2j) . which results in the presence of big 
gaps between these points and, hence, distorted values of the sphne interpolant. 

In Fig. El we show several plots of log^Q of the left-hand side of ()13p versus log^g t. The six 
lines in each group of lines of similar slope correspond to the six different choices fll6p - (PT|) 
of functions V, and the lines in each group come from the same value of 77 in i]Vd\i . Each 
on the "lines" in the figure in fact consists of 400 points (visually indistinguishable). The 
computer time spent on the CLP analysis is of the order of one minute per point (we used 
2^° Fourier coefficients to compute each of these points). 

We computed the regularity by performing linear regression on the points on graphs like 
the one in Fig. El in the regions where the points follow more or less a straight line. As one 
can see from this figure, for t close to 1 (i.e., log^^Q t ~ 0), the graphs for different functions are 
not straight lines, then as t decreases they form more or less straight lines, and as t decreases 
further, these lines level out. This behavior can be understood intuitively from (jl4j) - for 
t ^ 1 the high-fc Fourier coefficients are strongly suppressed by the factor (— 27r|/c|)'' e"^'^*''^', 
so the CLP method still does not "feel" the asymptotic self-similarity of the functions at 
small length-scales; in the other extreme, the leveling-out of the lines for very small t comes 
from the fact that in our computations we use a finite - albeit very large - number of Fourier 
coefficients. 

Remark 4.3. The "straight lines" in Fig. El are not really straight (which has been noticed 
in different contexts in jHElj)- We show this effect in Fig. ^1 which was created as follows. 
We took the six lines for 77 = 2 from Fig. El and for each of them we computed the slope 
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of the line as a function of the horizontal coordinate in the figure, log^^Q t. To compute this 
slope, we took each pair of adjacent points on the line and found the slope of the straight line 
connecting these points. The distance between two consecutive peaks in Fig. El is | log^Qcrd; 
more interestingly, as log^^gt becomes more negative, the lines tend to the same wavy line, 
until all lines reach saturation around log^o'^ ~ —4.5. 

4.4 Global Holder regularities — numerical results 

Table ^ summarizes our numerical results. The first column gives the map V used in the 
numerical computations (for the six functions V given by (jl6|) - (PT|) ). In the other columns 
we give the the values of the (global) Holder exponent k of the function R (representing 
the invariant circle as a graph in the (^, r)-plane), the advance map g, the conjugacy h and 
its inverse, h~^, coming from the (dynamics on) the golden critical invariant circle of the 
corresponding area-preserving twist map F. The notations used are the following: 1.85(15) 
stands for 1.85 ± 0.15, and 0.726(3) for 0.726 ± 0.003. Note that within the numerical error, 
k{R) = K{g), as expected. 

We also computed the Holder regularities of all big conjugacies H between each of the 
six functions hi, . . ., (coming from Vi, . . ., Vq) with all other h/s. We applied the CLP 
method to find that the regularity of all thirty functions H studied is 

/t(//) = 1.80 ±0.15 . (22) 

5 Discussion and conclusion 

In Sections 15.11 and 15.21 we point out some relationships between our results and previous 
studies related to universal scaling factors and singular measures. In the final Section 15. 3| 
we recapitulate our findings. 

5.1 Holder regularity and scaling factors 

Here we will explain how the scaling of the distances of closest returns of the iterates of a 
point gives bounds on the Holder regularity of some of the functions we study. Our analysis 
here is reminiscent of the analysis in jH Sec. 8.2]. 

We start by recalling the crucial observation of Kadanoff and Shenker Pfl |2] (see also 
jl5| Sec. 4.4]) of the existence of universal scalings in the distribution of the iterates of the 
Taylor- Chirikov map on the critical invariant circle 7 in neighborhoods of certain points 
of 7. Let ^rar ^ T staud for the value around which the iterates of the function g are most 
rarefied (in our notations 6'rar = |, while in j2] Orar = 0). Let ^^dcn ^ T stand for the value 
around which the iterates of the function g are most dense (in our notations ^rar = 0, while 
in j2j it is ^rar = |)- Since by Theorem 12.21 the function R is Lipschitz, around the points 
(^rar, -R(^rar)) and (6'den 5 -R(^den) ) 5 the iterates of any point on 7 under F are most rarefied, 
resp. dense. Shenker and Kadanoff found that the critical invariant circle in a neighborhood 
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of ^rar IS asymptotically invariant under simultaneous scalings in both 6- and r-directions, 
with scahng factors 

ao ~ -1.414836 (in 6) , /?o ~ -3.0668882 (in r) 

(see also the bounds on these values in Stirnemann [Sni)- This implies that, for large n, 



Po' ■ (23) 



The scaling around O^^n is a bit more complicated - it is called "step-3" scaling for obvious 
reasons: 

^7^"+^(^dcn) - ^dcn _l i?(^7'5"+^(^dcn)) - i?(^dcn) 



^?«"(^den)-^den " ' ' /2(^?«" (^den) ) - i?(^den) 

where the "step-3" scaling factors are 



/?3"' , (24) 



ttg ^ -4.84581 (in 9) , (3^^ -16.8597 (in r) . 

To understand heuristically why these scalings give restrictions on the Holder regularity 
of i?, set := g'^"'+^ (O^ar) — ^rar, := R{g'^"+^{Oj.a.r)) — -R(6'rar) for somc large value of n. 
Then if the local Holder exponent of i? at ^ = 6'rar is n, we will have |Ar| ~ |A^|'^. If 
the graph of R is asymptotically invariant around (6'i.ar, R{Ora.r)) with respect to the scalings 
we will have |/5o Ar| ~ lao^^^. "Dividing out" the last two relationships, we obtain 



\I3q\ ~ |ao|'^, i.e., k, ~ . This argument (which can easily be made rigorous) implies 

that the (global) Holder exponent of R does not exceed ^ 3.22945. The scaling ^ 

yields a tighter bound on the Holder regularity of R: 

^i{R) < ^^^^ ^ 1.7901 . (25) 



Note that the fact that the scaling ()24|) is "step-3" (as opposed to "step-1") is irrelevant for 
the bounds on the Holder regularity. 

To obtain bounds on K{h) and n{h^^), we use Lemma 8.1 from j4j, which says that if the 
function h conjugates fi and f2, ho fi = f2oh, and if for some sequence of positive integers 
Qn the functions fj (j = 1, 2) behave in a neighborhood of the fixed point 6'fix = /i(6'fix) of h 
as follows: 

/f"(M=^fix + C,r;-- + o(r77") 

for some constants rjj and Cj, then K(h) < Applying this to the definition of h and 

using the well-known fact that {Qn ctg) ^od 1 < Ca^, we obtain the bounds 

^(h) < If^sK^ _ 0.721125 , ^i(h-^) < i^^^^ ^ 0.91478 . (26) 
" fog|o-G| ^ ' - log lag 'I ^ ^ 

A comparison with Table [T] suggests that these bounds are saturated. 
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5.2 Conjugacies and singular measures 

The functions whose Holder regularity we study are defined through high iterates of maps. 
For example, the graph of the function R defined by (0)) is nothing but the the critical invari- 
ant circle 7 of F which is is filled densely by the iterates F"-[6o, tq) of some point {60, tq) G 7. 
Here we discuss how some characterizations of the singularities in the distribution of the 
iterates of F on 7 are related to the Holder regularity of some of the functions considered. 

Hentschel and Procaccia |27j pointed out the importance of the generalized (Renyi) di- 
mensions D{q) of a singular measure for dynamical systems; these quantities have been 
defined previously in the context of probability theory by Renyi Halsey et al in their 
seminal paper jS] related heuristically the Renyi dimension of a singular measure to the 
spectrum of singularities f{a). We recall that f{a) is the Hausdorff dimension of the set Ea 
of points where the measure has singularity of strength a. The spectrum /(a) is a function 
supported on the interval [amin^ttmax], where = D{oo), resp. ctmax = D{—oo), describe 
the scaling behavior of the measure in the region where the measure is most dense, resp. 
most rarefied. 

Let (^0) '"o) be an arbitrary point on the critical invariant circle 7 of the area-preserving 
twist map F. Then the distribution of the iterates in a very long orbit, {F"'(6'o, ro)}^0' 
approaches as ^ 00 the "density" of the measure on 7 that is invariant with respect 
of the restriction of the map F onto 7. (We put "density" is quotation marks because for 
singular measures this is not a function, but a set of Dirac (5-distributions.) This invariant 
measure on 7 induces an invariant measure /ig of the map g on T. It is easy to see that (|7j) 
implies that 

h-\e) = / d^g 

Jo 

(for an appropriately chosen ( in the redefinition of h as in Sec. 12. This relationship 
imphes that the spectrum of singularities /(a) of the measure /ig is the same as the Holder 
spectrum /h(«) of the function h~^. By definition, /h(«) is the Hausdorff dimension of the 
set where the local Holder exponent of the function is equal to a; for a readable account we 
refer the reader to Jaffard [HH]- The (global) Holder regularity /t(0) of a function (p is equal 
to the lowest end, amin, of the support of the Holder spectrum, /h(«), of 0. 

Osbaldestin and Sarkis jlHl applied the method of to determine numerically the func- 
tions /(a) and D{q) of the invariant measure fig coming from the distribution of iterates of 
the Taylor- Chirikov map F on the golden invariant circle. They found that 

ttmin = D{oo) ^ 0.915 , Omax = D{-0o) ^ 1.387 ^ . 

Comparing with the values in Table d the reader should recognize that their a^i^ is nothing 
but our K{h~^), while amin is equal to the inverse of the regularity of the conjugacy h. 

Buric et al |1TI studied numerically the Taylor-Chirikov map and the map (P) with 
V{6) = I sin27r^ -|- | sin 4-7?^, for rotation numbers with continued fraction expansions of the 
form [S, 1°°] := [S, 1, 1,1,...], [S, 2°°], [S, 3°°], [S, 4°°], where S stands for some short string 
of positive integers. They found that /(a) and D{q) depend only on the tail but do not 
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depend on the initial part S as well as on whether the Taylor- Chirikov map or the other 
map was used in their numerics. 

Other papers related to numerical computations of singular measures on critical invariant 
circles of area-preserving twist maps are Shi and Hu |13I33], where the methods of |3j were 
used, and Hunt et al where the authors used the thermodynamic formalism developed 
in jlHl to compute the information dimension -D(l) of the standard map for different rotation 
numbers. 

5.3 Conclusion 

We computed accurately the golden critical invariant circles for six twist maps of the form 
and the global Holder regularity k of some functions related to the dynamics on these circles. 
Our numerical experiments lend credibility to Conjectures 12.31 12.41 and 12.51 concerning the 
universality of the regularities of the functions R, g, h, and H (see Tabled] and (I221))- 
Yamaguchi and Tanikawa |17] found numerically that the golden invariant circle (given by the 
function R) of the Taylor- Chirikov map is differentiable but R' is not of bounded variation; 
our studies significantly narrow the numerical bounds on k{R) for this and for other maps. 

Our results seem to indicate that the regularities of R, h, and saturate the upper 
bounds and (PBj) coming from previous studies of scaling exponents. 

Our finding that k,{H) is greater than K{h) and n{h~^) by a comfortable margin (cf. Con- 
jecture 12. 6p has an interesting consequence. As discussed in Sec. 15.21 the Holder regularity 
of h and is different at different points, and for each a G (amin, «max), the set (where 
the pointwise Holder exponent of is a) has Hausdorff dimension fH{<^) strictly between 
and 1. Previous numerical studies indicated that fuiot) are the same for different maps F. 
Our finding shows that the "irregularities" of functions h coming from different maps F are 
interspersed in the same way in [0, 1] for all twist maps studied. Note that this does not 
mean that for a certain value of a the sets E^ are the same for different F in the same 
universality class - only the way all sets E^ for different a are interwoven is universal. 

It would be interesting to apply wavelet-maxima methods for pointwise regularity [^ HU] 
(see also the rigorous analysis in ISHl) to the problem studied in this paper and to compare 
the results of the wavelet analysis with the results about the singular invariant measures. 

As a by-product of our studies, we have computed millions of Fourier coefficients of the 
functions /i, and noticed some self-similarity properties that to the best of our knowledge 
have not been observed before. Presently we are working on understanding these properties. 
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Tables 



F with: 


k{R) 


1^(9) 


K{h) 






1.83(9) 


1.83(9) 


0.722(1) 


0.92(1) 


V2 


1.79(6) 


1.75(9) 


0.721(1) 


0.92(1) 




1.83(4) 


1.84(3) 


0.724(2) 


0.93(2) 




1.86(8) 


1.86(8) 


0.722(1) 


0.92(1) 




1.85(5) 


1.85(5) 


0.724(2) 


0.93(1) 




1.85(15) 


1.88(12) 


0.726(3) 


0.93(2) 



Table 1: Regularities of the functions R, g, h, and h ^ for the golden critical invariant circles 
of different maps F. 
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Figure captions 



Caption to Figure [TJ 

Critical invariant circles, r = R{6), of the maps corresponding to the maps Vi, V2, . . ., 
Vq given by (fTm) - (PT|) (Vi = thin solid line, V2 = thick solid line, V3 = dotted line, V4 
= thin dashed line, V5 = thick dashed line, Vq = dot-dashed line). 

Caption to Figure [2J 

"Periodized" advance maps g — Id (notation same as in Fig. 
Caption to Figure |Ht 

"Periodized" conjugacies h — Id (notation same as in Fig. [T}. 
Caption to Figure |3t 

"Periodized" inverse conjugacies — Id (notation same as in Fig. Q). 
Caption to Figure 

Zooming in the graph of the function h — Id corresponding to the map V2 p7|) . 

Caption to Figure |Ht 

"Periodized" big conjugacies H — Id. 



Caption to Figure [71 

Plot of logio (h^^ 



versus log^g k, where h corresponds to the map F coming form 



the function V3 (fH^ll . 
Caption to Figure |St 



Plot of \ogiQ [g — Id 



and log 



10 



/i-i - Id 



versus log^^Q k, for the same map F as 



in Fig. [7| The impulses correspond to {g — Id), and the dots above them to {h ^ — Id). 



Caption to Figure IHt 

Plots of log 



10 



9 ^-ts 



(I) 



L°°(T) 



versus log^^Q t for the functions K = {h — Id) for 



the twist maps coming from Vi, . . ., Vg, for = 2 (shallowest lines), rj = 3, and 77 = 4 
(steepest lines). 

Caption to Figure llOt 

Slope of the lines on Fig. El as a function of log^Qt (see the text). The notation is the 
same as in Fig. [H 
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Figure 2: 
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Figure 3: 
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Figure 4: 
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Figure 10: 
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